Free-energy transition in a gas of non-interacting nonlinear wave-particles 
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We investigate the dynamics of a gas of non-interacting particle-like soliton waves, demonstrating 
that phase transitions originate from their collective behavior. This is predicted by solving exactly 
the nonlinear equations and by employing methods of the statistical mechanics of chaos. In par- 
ticular, we show that a suitable free energy undergoes a metamorphosis as the input excitation 
is increased, thereby developing a first order phase transition whose measurable manifestation is 
the formation of shock waves. This demonstrates that even the simplest phase-space dynamics, 
involving independent (uncoupled) degrees of freedom, can sustain critical phenomena. 

PACS numbers: 05.00.00,47.35. Jk,02. 30. Ik 



The methods of statistical mechanics have permeated 
physics as a whole including modern areas of determinis- 
tic chaos [l[ , complexity Q , and nonlinear waves 0, 0> S] • 
Among its basic notions, it is well known that the ther- 
modynamics associated to a gas of non-interacting par- 
ticles is trivial, being characterized by a constant free- 
energy with respect to control parameters. Conversely, 
in order to have cooperative phenomena, such as phase- 
transitions or more complex processes, the intervention 
of some sort of interaction potential must be called for. 

Within the context of nonlinear wave propagation, a 
challenge to the generality of this picture could come 
from solitons, i.e. exact solutions of integrable nonlin- 
ear partial differential equations (PDE), which display 
well-known particle-like behavior. Stemming from the 
numerical experiments of Fermi-Pasta-Ulam (FPU) on 
the equipartition of energy in nonlinear chains [§], soli- 
tons have found applications in area as different as Bose- 
Einstein condensation (BEC) [7], nonlinear optics @, 
fluid- dynamics (Toj , solid state physics [IH , general rela- 
tivity [12j | and many others. While relying on large (non- 
perturbative) nonlinear effects, solitons conserve their 
number and spectral parameters (eigenvalues) evolving 
without any practical interaction, albeit some displace- 
ments in their collisions. The isospectrality allows for 
reducing the infinito-dimensional phase space associated 
to the global wave-function to a simple one where N in- 
dependent degrees of freedom corresponding to N soli- 
ton particles are effective. On this basis one can argue 
whether the statistical description of an ensemble of non- 
interacting solitons is actually as trivial (i.e., with a con- 
stant free energy and no critical phenomena) as inferred 
from the statistical mechanics of free particles, or else ex- 
hibits cooperative phenomena as the number of particles 
(degrees of freedom) grows. In this Letter, we address 
this issue by investigating the statistical mechanics of 
a gas of soliton particles. By exploiting ideas from the 
thermostatistics of chaos [lj , we demonstrate that a suit- 



ably defined free-energy undergoes a metamorphosis as 
the number of solitons grows, changing from a constant 
to a function that supports a first-order phase transi- 
tion. Correspondingly, the overall wave-function devel- 
ops a steep front (shock), which results from the scaling 
properties of the soliton velocity distribution, in turn de- 
termined by the input excitation. Hence the shock for- 
mation can be interpreted as a cooperative process of 
several solitons. Therefore even the simplest phase-space 
dynamics, where particle-like waves evolve independently 
from the others, could result into phase-transitions. 

We make reference to a universal integrable model, 
namely the defocusing nonlinear Schrodinger equation: 

and by means of the inverse scattering transform (1ST) 
[ill, we calculate the field evolution after the so-called 
dark initial condition ip(x,Q) = — ptanhx, with p inte- 
ger. The latter is experimentally accessible andpartic- 
ularly important in BEC and nonlinear optics [J S, @] , 
where the occurrence of dispersive shock waves has been 
recently observed and discussed [li , H, 3, 17, 3, [3] , 
though in regimes not involving dark disturbances and 
the soliton-driven wave-breaking discussed here. The 
constraint p=integer, corresponding to a reflectionless 
potential, allows us to investigate a novel scenario where 
the dispersive shock dynamics is determined solely by 
solitons that are embedded in the input, with radiation 
waves playing absolutely no role. Unlike previous ap- 
proaches to shock waves (hydrodynamic limit, Whitham 
averaging) involving different degrees of approximation 
HI 0, El HI HI, we derive exact solutions of Eq. Q, 
which enable us to contrast the smooth dynamics deter- 
mined by a small number N of solitons with the strongly 
nonlinear case (p, N ^> 1). In the latter regime, this 
approach allows us to introduce a free energy that re- 
flects the spectral distribution of solitons in the (invari- 
ant) eigenvalue space, using a non-canonical measure 0]. 
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In this respect, our approach shows no similarities with 
the Gibbsian statistical mechanics of soliton-bearing sys- 
tems where nonlinear excitations are thermally controlled 
[H> 0] j or with the kinetic theory of random distributions 
of colliding solitons Our formulation is general and 
applies to every integrable system, such as integrable ver- 
sions of FPU. 

Spectral transform. - The solution of Eq. ([I]) is de- 
rived by means of 1ST in two steps: (i) by calculating 
the spectral transform of the system; (ii) by solving the 
inverse scattering problem. It is worth to remark that 
the solution of Eq. |T]) with initial value tpo = tp(x,0) 
represents a far more general problem than the calcula- 
tion of the generic N— soliton solution of Eq. (1). The 
former, in fact, calls for the evaluation of the spectral 
transform S{ipo} of tpo and the solution of the inverse 
problem for S{tp }, while the latter needs just the cal- 
culation of the inverse problem for an unknown tp whose 
spectrum contains only solitons. In order to highlight the 
main physical results, we defer unessential mathematical 
details to an extended publication, while referring to Ref. 
[p| for the 1ST method. The spectral transform S{tpo} 
is defined as: 



S{tpo} = {— oo < A < oo. 
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(2) 

with the reflection coefficient b/a — R, A„ correspond- 
ing to simple zeros of a (for Im(A) > 0) and c„ = 
&(A„)/a(A„), being A a spectral parameter and a(A), &(A) 
transition coefficients relating the Jost solutions of the 
direct scattering problem . We calculate a and b cor- 
responding to tpo, obtaining: 



b= fc-|r(z|)| 2 



■A.r(i + p)r(A- P )' 2r(i-p)r( P )' (3) 

When p is integer, R = and tpo is a reflectionless 
potential containing N p = N(p) = 2p — 1 solitons. The 
spectral transform S{tpo} then reads: 



S{M = W) = 0; A±„ = ±2^p 2 -{p 



C±n = (-1) 



j^n+l ]lj( S ±n, 



Uj=i±J S ±n - Sj) 



,n = 0, 



. M}, (4) 



with s n = A„ + 2i(p - n), j G [-M, M] and M = p — 1. 

Inverse scattering problem. - The ^-evolution of the 
wave function tp is found from the solution of the inverse 
problem [131 ] and, for the spectral transform represented 
by Eq. ((3]), reads: 



ip(x,z) = p 
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where di — i^/a/sie 
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-M, M], Vj = 2(p — j) and A is the matrix 
A, r 1 -■ HV 3 



Sji + i- 



(6) 



0.33 
P0.16 


0.5 j 

0.4 
0.3 

z 

0.2 
0.1 

o! 



f a } 1 

xo 
-1 



(b) 



-0.5 0_ 0.5 

A/V3p 



1 



w 



-2.5 



2.5 




FIG. 1: (Color online) 3-dark soliton solution corresponding 
to tpo — 2tanhj:: (a) Eigenvalue measure P (pj/e for e — *• 0); 
(b) numerically calculated free energy X vs. j3; (c) Pseudo- 
color level plot of intensity IV'I 2 ; (d)-(e) Profiles of frequency 
S x and intensity \tp\ 2 at different z . 



Equations (gj) and ([5]) describe the soliton evolution cor- 
responding to the input tp . Note that A„ represent the 
invariant velocity of the n-th soliton. 

Soliton-gas thermodynamics. - To derive a soliton- 
gas thermodynamics, we begin by decomposing the spec- 
tral transform S{ipo} as the sum of spectral transforms 
of individual solitons (normal modes) : 
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with \j and Cj given by (|3]). A remarkable property of 
integrable PDE is that each conserved quantity of the 
system, such as the energy or the momentum, follows 
the decomposition ([7|) and can be expressed as the sum of 
conserved quantities of individual solitons. By exploiting 
this property, we evaluate the Hamiltonian H: 



H 



dx 



dtp 

dx 
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(8) 



whose explicit summation yields H — ^p 2 (p 2 
the Hamiltonian of tp — tpo = — ptanha;. Equation (jSJ) 
states that the system can be regarded as a gas of non- 
interacting particles possessing a purely kinetic Hamilto- 
nian of the form H = Y)j mVj/2, with mass m = 2/3 

3/2 

and equivalent velocities Vj = ty . The latter governs 
the observable dynamics in the "dual space" of the wave- 
function tp, through Eqs. ([5][6]). 

The thermodynamics of such soliton gas can be studied 
by means of escort distributions [l[. Given an arbitrary 
measure, whose representative microstates have probabil- 
ity pi, escort distributions are constructed from the Lie 
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transformation group Pf. = f (pk ,13) = pf / Ylj Pj ■ The 
connection with statistical mechanics is derived by defin- 
ing the partition function Z(J3) = J^jPj — ex P(~ 
correspondingly Pj = cxp(4' — (3£i), £i = — lnp, being 
the energy of the i— th microstate , and /3 = 1/T > 
the inverse "temperature". The Helmholtz free energy 
is given by F{j3) — ^/(3. We define the phase-space as 
the range [—1,1] of A/2p. As a principal measure of the 
soliton gas, we consider its eigenvalue distribution \i/2p, 
which is at the basis of the evolution in the dual wave- 
function space through Eq. ([5]). We then partition the 
phase-space in N e = 2/e boxes of size e, and define pi 
as the normalized Q^yPj = 1) probability of finding an 
eigenvalue in the interval [—1 + (i — l)e, — 1 + ie] with 
i = 1,2, >.„N e . Following the thermostatistics of mul- 
tifractals [l[, we study the thermodynamic limit of P, 
defined by 




(3F=X 



lim 



y 



(9) 



being V = — In e the "volume" . Correspondingly, when 
e — * 0, the partition function scales as Z ~ e x . 

Thermodynamic limit for small particle densities. - 
We begin by calculating A" when the amplitude p is small 
enough so that the N p eigenvalues Xj are well separated 
in the spectrum. This situation corresponds to a low 
density soliton gas. In this case, the probability pj of 
finding an eigenvalue in the j— th box of width e is either 
or 1/Np. Therefore, in the limit e — > 0, the non- vanishing 
Pj scales as e , since each eigenvalue occupies a single 
point in the gap |A/2p| < 1. 
reads 



FIG. 2: (Color online). (a) Eigenvalue measure P (for e — > 0) 
versus A/2p for a dense gas; (b) Free energy X versus inverse 
temperature ft. 



same position x = 0) to the asymptotic stage where they 
are well separated. As stated above, this relies on the 
absence of phase transitions in soliton gas (Fig. QJi-b), 
owing to a constant free energy landscape T = 0. In 
summary, cooperative phenomena and shocks are pro- 
hibited at small input amplitudes p, or equivalently, for 
low density soliton gases where T = 0. 

Thermodynamics for dense soliton gases. When 
the number of eigenvalues A grows large for p ^> 1, the 
soliton ensemble turns into a dense gas and the eigen- 



value distribution X/2p — ±y 1 — (1 — n/p) 2 becomes a 
continuous function of the variable n/p. The density of 
the positive eigenvalues T>(X/2p) is readily seen to be 

T>(y) 
that 



yj2\J\ — y 2 and, correspondingly, it is found 



Vi = v/i - (J - i)V - V^J 2 



(12) 



The partition function then ^or e — > 0, we identify three different scaling regimes in 



Eq. (fT2~|) : (i) pj cx e for A = 0; (ii) pj cx e for A = p; 
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(iii) pj cx ^/e for A 
scaling as 



2p, with the partition function Z 



(10) 

This yields a very simple thermodynamics, where both X 
and J 7 are zero, with the distribution function P of the 
eigenvalues, determined as the limit of Pj/e for e — * 0, 
being a set of Dirac delta (Fig. [T^,b). The resulting 
wave evolution is smooth, without the formation of any 
singular behavior. We specifically address the case of p = 
2, where the results of the 1ST analysis are manageable to 
be reported in simple closed form. Equation (|5|) yields the 
particular 3-dark soliton solution relative to eigenvalues 
0, ±2^3 as a function of X = 2x, ( = 4V3z G [0, oo): 



Z 



(N e - 3)e?+2e^ 2 + 



(13) 



with two boxes near A = ±2p scaling like e' 3 / 2 , one near 
A = scaling like e 2 ^ and the remaining N e — 3 scaling 



like . In the limit 



0, we can actually neglect the 



contribution of e 2/3 (as the probability pj = for A = 0). 
The partition function then reads Z ~ e 



This implies X — min[/3 — 1, (3/2] and leads to: 



X = 



4(coshC + cosh A + «\/3sinhC) sinhA 




f3<(3c 
(3>(3 C 



(14) 



(11) 



— 3 — 4 cosh £ cosh A — cosh 2 A 

Figures [TJ:-e display the level plot evolution of the inten- 
sity or density (in BEC) |i/>| 2 , and snapshots of |^| 2 and 
S x = dS/dx = lm{ip x /ip}. The generation of three dark 
solitons by ipo does not lead to the formation of steep 
fronts during propagation: each soliton slowly splits up 
from the others throughout the process, from the ini- 
tial overlap state (i.e., at z = the particles occupy the 



being (3 C = 2 the critical point such that (3 C — 1 = (3 c /2. 
Equation (TTJJ states that the system undergoes a first- 
order phase transition as the inverse temperature f3 is 
varied, since the free energy (3T = X is continuous but 
not differentiable at (3 C . As seen in Fig. [2p, X and P cal- 
culated numerically after the expression of the eigenvalue 
X n support this picture. Such a phase transition is the 
main mechanism leading to shock formation in the dual 
space of wave dynamics. In fact, owing to the different 
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formation of a dispersive shock wave involves only soli- 
tons. In spite of the fact that solitons behave as non- 
interacting particles with well-defined parameters, they 
do not follow the canonical behavior of systems with non- 
interacting degrees of freedom. Rather their cooperation 
leads to measurable critical phenomena characterized by 
a free energy that develops a singularity when the num- 
ber of solitons grows sufficiently large. This finding can 
stimulate both novel experiments and new ideas in the 
field of statistical mechanics of nonlinear waves. 



FIG. 3: (Color online), (a) Level plot of intensity IV"! 2 , and 
(b) snapshots of frequency S x and (c) intensity \ip\ 2 , for input 
ipo — 30tanhir. 



scaling of pj , wave-particles with close but opposite soli- 
ton velocities Xj tend to be accumulated at gap edges; 
starting from the initial overlap state, such process gen- 
erates a nearly discontinuous variation in the overall fre- 
quency S x - To illustrate this dynamics, we calculate the 
evolution of a dense soliton gas corresponding to p = 30 
(Fig. [3]). A dramatic variation of S x is observed for 
z = z* ss 0.013 (Fig. [3p), with a shock (steep front) 
in x = eventually originating from a compressional-like 
wave (S x > for x < and viceversa), beyond which fast 
oscillations in x appear. Each oscillation is indeed asso- 
ciated with a wave particle (dark soliton), with solitons 
splitting up in pairs with opposite velocities around a 
central black (A = 0) soliton (Fig. [3^). Correspondingly, 
the intensity I'i/'l 2 becomes singular in z* (Fig. [2b). The 
dynamics of ip in the dense case (p = 30) is dramatically 
different if compared to the 3-soliton case (see Figs. Hh-d 
and[3J>c): the absence of cooperative phenomena when 
the number of soliton is small, in fact, does not lead to 
any singular behavior (shock wave) in the field propaga- 
tion along z. We emphasize that, despite solitons tend 
to split up from the beginning (z = 0) due to their dif- 
ferent velocities, their collective behavior results, owing 
to the compressional wave feature, into a dramatic focus- 
ing of the dark notch intensity, which turns out to be an 
easily measurable signature of the critical behavior. Re- 
markably this occurs while: (i) no radiation is involved; 
(ii) the phase space supporting the gas dynamics is com- 
posed by non-interacting particles. In summary, a dense 
gas of noninteracting soliton wave-particles turns out to 
support critical phenomena. The latter originates from 
a singularity of the free-energy, which undergoes a meta- 
morphosis as the number of soliton grows and develops a 
first-order phase transition. As shown, the relevant order 
parameter for the soliton gas transition is the shape of the 
function /3J 7 , which changes with the input excitation. 
In conclusion, we have illustrated a scenario where the 
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